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We discuss a hopping model of electrons between idealized molecular sites with local orbital 
degeneracy and dynamical Jahn- Teller effect, for crystal field environments of sufficiently high sym- 
metry. For the Mott-insulating case (one electron per site and large Coulomb repulsions), in the 
simplest two-fold degenerate situation, we are led to consider a particular exchange hamiltonian, de- 
scribing two isotropic spin-f/2 Heisenberg problems coupled by a quartic term on equivalent bonds. 
This twin-exchange hamiltonian applies to a physical regime in which the inter-orbital singlet is the 
lowest-energy intermediate state available for hopping. This regime is favored by a relatively strong 
electron- phonon coupling. Using variational arguments, a large- n limit, and exact diagonalization 
data, we find that the ground state, in the one dimensional case, is a solid valence bond state. The 
situation in the two dimensional case is less clear. Finally, the behavior of the system upon hole 
doping is studied in one dimension. 
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I. INTRODUCTION 



Ions or molecules with degenerate electronic orbitals coupled to local vibrations can display a rich variety of 
phenomena. |5J The most notable of these phenomena is the so-called static Jahn- Teller (JT) effect, i.e., the system 
gains energy by lowering the symmetry through a static distortion, thereby removing the electronic degeneracy. For 
finite systems, like a molecule, quantum mechanics can play an important role in that, if the frequency of the relevant 
vibronic modes is high enough, it can dynamically restore the original symmetry through a quantum superposition of 
"degenerate" distorted configurations (dynamical JT effect). 

Many systems containing JT ions or molecules are found to be insulators with quite a rich variety of spin and orbital 
ordering effects. Such a richness of exchange effects is due to the interplay between the spin and orbital degrees of 
freedom. (!]] 

The simplest cases can be modeled by the so-called E <g) e JT-coupling, i.e., at each site r a two-fold degenerate 
electronic orbital is coupled to a two-fold degenerate local vibron. Q Although conceived as the simplest workable 
model, such a situation is not totally unrealistic, as higher molecular degeneracies can often be reduced, by the crystal 
field of the solid state environment, to residual doublets. (If the crystal field symmetry is very low, the degeneracy 
will be completely removed). 

If exactly one electron sits on each JT site, and if the Coulomb repulsion is very large, the system will have 
an insulating gap, due to Mott-Hubbard correlations, and the low-energy spectrum will consist of spin and orbital 
excitations. For the E e case, which we are going to consider, the physics of the insulating state is governed by 
an exchange hamiltonian involving, besides the actual spin S r , also a pseudo-spin- 1/2 variable T r , representing the 
orbital degrees of freedom. Unlike the case of a non-degenerate orbital, where antiferromagnetism is favored, one 
finds here a rich phase diagram, depending on the energy of the two-particle intermediate states available for virtual 
hopping. H If Hund's rule coupling is the dominant effect, the lowest doubly occupied site is a spin triplet, and spin 
ferromagnetism turns out to be favored. Jj[] Many experimental systems, jij including possibly the case of TDAE-C@o, 
P,PJ may fit into such a scenario. In the present paper we are going to consider a different regime, i.e., when the 
coupling to vibrons is strong enough that an inter-orbital singlet state becomes the lowest-energy doubly-occupied 
virtual state. We are going to argue that new interesting phases can show up in this case. The exchange hamiltonian 
for this case can be shown to be related to the following spin-model: JTo| 

H ST = -J ]T (2S r • S r , - ~)(2T r • TV - ~) . (1) 

<rr'> 

We will show that, in this case, valence bond (VB) states tend to be favored, at least in low enough spatial dimensions. 

The paper is organized as follows. Section [n| introduces the hopping model that we will consider. In section III 
we will discuss in detail the hamiltonian -ffsT m Eq. ( pi). In section ^ we study the behavior of the model upon 
introduction of holes, in one dimension. Finally, sectionM contains our conclusions, and some open issues. 



II. THE MODEL 



Suppose that at each site r of the lattice we have two degenerate electronic orbitals and two degenerate phonon 
modes, of frequency ujq, both labeled by an orbital quantum number I — ±. Let cL(r) be the creation operator for 
an electron with spin a and orbital quantum number I at site r, and b, (r) the phonon operator. In a standard E <g> e 
Jahn- Teller problem, the on site linear electron-phonon interaction term can be written as Q 

ffe-ph(r) = uj J2i b l b l + V 2 ] + W2){T+ [bl + b + ] + H.c.} , (2) 

;=± 

where T = (1/2) ^ Q J^i v c i a {&)ll' c l',at is the orbital (pseudo) spin operator, <x are Pauli matrices and we have omitted 
the obvious site-label. We will also consider electron-electron terms, which include a Hubbard intra-orbital repulsion 
U, an inter-orbital repulsion Ui, and an exchange (Hund's) coupling Jh, 

H e -e(r) = U } ni^n^i + Uj }^ n + ^n-^> - J«S + • S_ , (3) 

l— ± cr.er' 

where S { = (1/2) £ , cj ja) ci a i is the spin operator for the ^-orbital. 
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Assume now that there are (small) electronic hopping processes between the different sites described by a standard 
tight-binding hamiltonian. The total hamiltonian is therefore written as: 

H = -t E[ c tW^(/)+H.c] + £iT site (r), (4) 

<rr'> I, a r 

where we have neglected overlap integrals between orbitals with different orbital quantum numbers, and 

H sitc (r) = fle_ ph (r) + flo-e(r) • (5) 

We will be interested in the limit in which double occupancy of a site is strongly inhibited by the repulsive on-site 
interactions. The low-energy physics is then described by an effective hamiltonian in which the double occupancy 
is only treated as a virtual process to second-order in the hopping matrix element t. For a two-fold degenerate 
orbital, there are six possible two-particle states to be considered as intermediate states: an inter-orbital singlet 
1 5) = \S = 0,T Z = 0) = 2- 1 /2( c ^ T c^ i _ c t_ )i c^ T )|o), two intra-orbital singlets |±) = \S = 0, T z = ±1) = 4 T 4 j0), 
and an inter-orbital triplet \T) = \S = 1,T = 0). For the case of one electron per site, n = 1, the resulting effective 
hamiltonian involves two spin-1/2 operators at each site, the spin S and the orbital pseudo-spin T. The derivation is 
standard, [|||||| yielding 

ffeff = 2 tt Js ~ Jt + 2J± ) S '' ' Sr ' + ( 3Jt " Js ^ Tr ■ T '" + 2 ( Js " J ±) T r T r'] 

[4(Js + Jr)(S r ■ S P /)(T r • TV) + 8(J± - J s )(S r ■ S r ,)(T r z T 7 z ,)} . (6) 



<rr'> 



The couplings Js, Jt, and J±, generalization of the famous antiferromagnetic coupling At 2 /U of the non degenerate 
Hubbard model, are essentially given by t 2 a j AE a where AE a is the energy difference between a doubly occupied 
site of type a = S,T,± and two singly occupied sites, and t a is an effective hopping matrix element reduced by 
electron-phonon interactions (Ham's factor). || 

In absence of electron-phonon coupling (g = 0), t a = t and AE^ = Ui + (3/A)J H , AE^ — Ui — (1/4)J#, 

AE^ = U. ||§ The triplet state is therefore the lowest energy virtual state (under the realistic assumption Uj < U), 
and Jt the largest coupling in the problem. Spin ferromagnetism is the natural outcome of such a scenario. |l],|3| The 
Jahn- Teller coupling can change this picture considerably. First of all, standard polaronic effects lead to a decrease 
of the energy of the singlet states, primarily the inter-orbital singlet \S), with respect to the triplet states \T). [f|@||] 
Second, the effective hopping matrix clement t a turns out to be larger for the singlet \S) than for the other states. 
Ref. H gives the results for the J's obtained from perturbation theory in <?, accurate to order <? 4 : 

7 = t' 2 (l+g 4 /2) 

' 3 U* + (3/4) J H - (,g 2 /2 - 7. 9 4 /16H 1 ' 

Jt — 

J± = 



Ui - (1/4) J H + (,9 2 /2 - 
t' 2 



U - (3 ff 4 /i6Vo ' 

In the opposite strong electron-phonon coupling case, g — > oo, it was shown in Ref. [g| that, again, the most important 
doubly-occupied state is the inter-orbital singlet IS*), and the model can be mapped into a single-band hamiltonian 
with a spin-1 variable attached to each site. |l^] This also leads to a situation in which, effectively, Js is the only 
coupling left in the problem. |l(]] 

The mean field phase diagram of the model @ has been studied in Ref. || , together with a few special points where 
exact solutions can be obtained in D — 1. As it turns out, there is a region of the phase diagram, not considered so 
far, where, as we will argue below, valence bonds (VB) are favored. The physical regime in which this region is of 
relevance is precisely the JT limit in which Jg is the largest coupling in the problem, i.e., the two-particle inter-orbital 
singlet has lower energy than the corresponding triplet. 

More in detail, let us assume that Js > Jt and define J = Js — Jt- In a bipartite lattice the effective hamiltonian, 
after performing the canonical transformation — > (— l) r T^, can be written as: 

H eS = H ST + AH (8) 

where 
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Hst — J 2_^i I s - ■ S -' + T r ■ T r' ~ 4(S r ■ S r ,)(T r ■ T r ,)} , (9) 

<rr'> 

and 



AH = AJ T J2 [ T r T r' ~ ( 2S r ■ S r / + 1/2) ■ Tjb] 

<rr'> 

+4J± ( 2S *r • &r> - l/2)(T r z T* + 1/4) . (10) 



<rr' > 



In the next section we will concentrate on Hst, showing first of all that it certainly has a VB-solid ground state 
in D = 1, and suggesting that the same might be true in D = 2. The presence of a gap in the excitation spectrum 
implies the robustness of such a phase with respect to a small enough perturbation AH. 

Next, we will consider the question of the behavior of a valence bond state upon doping. The resulting generalization 



of the t — J model will be discussed, for the D = 1 case, in section IV 



III. ONE ELECTRON PER SITE: THE SPIN HAMILTONIAN. 

As previously discussed, when there is exactly one electron per site (n = 1) and one assumes the inter-orbital singlet 
\S) to be the only important intermediate two-particle state, the exchange hamiltonian governing the spin and orbital 
degrees of freedom can be written as 

H ST = - J ( 2S r ■ S r' ~)(2T r • T r , -\) = J £ (P™ l)(P rV ? - 1) , (11) 

<rr'> <rr'> 

where pj^) = 2S r -S r / + 1/2 is a permutator between the spins at sites r and r', and similarly p( T ) for the pseudo-spins. 
This form will be particularly useful in the following. 

This model shows an obvious spin and pseudospin rotation invariance described by the SU (2) £g) SU (2) group. 
However, the full symmetry group displayed by ([ll]) is much larger and includes SU(4). Indeed, one can show that, 
in a bipartite lattice, H ST commutes with the 15 operators S a = £ r S?, T a = £ r T r a , and = 2Y /r (-l) r SfT" 
which generate the Lie algebra su(4) pj] |. Due to the staggering in the L^ v operators, these generators do not 
commute with the spatial symmetries which interchange the two sublattices. Therefore, the full symmetry group is 
even larger than the aforementioned SU(4). Finally, we notice that this hamiltonian differs from the exactly soluble 

( 9^ (T} 

SU(i) Sutherland model H = JE<rr'> Kr'Prr' obtained by a natural generalization of the Hcisenberg hamiltonian. 
Rather, it is possible to show that the S'[/(n)-invariant model introduced by Affleck |l2) 

\ E W^(r'), (12) 

<rr' > 

is unitarily equivalent to Hst/^J for n = 4. Here Sp and S@ generate, respectively, the fundamental and anti- 
fundamental representations of SU(n). Q 

In D = 1, Affleck has shown that H n , in the limit n — * oo, has two simple VB-solid ground states, illustrated in Fig. 
0. Moreover, the n = 3 case has also been shown to have a dimerized ground state, although with a large spin-spin 
correlation length. |l3| As the tendency to dimerization is likely to increase with n, Jl2| we expect that the n = 4 
case, of interest to us, will also have a VB-solid type of ground state. This expectation will be explicitly shown to be 
correct in D = 1. 

Let us first show that singlet bonds between S and T variables are very natural objects to introduce in our problem, 
for any D. If we denote by (rr')s = 2~ 1 / 2 (trlr / — |rtr')s a singlet bond formed between the spins at r and r', it is 
easy to show that Q 

(P$ - l)(rr') s = -2{rr') s 

{P^-l){rr')s{ll')s = {r'l)s{rl') s . (13) 

Similar results apply to p( T ) — 1 when acting on valence bonds formed from the pseudo-spin variables, (rr')x- 
If we denote by (rr') = (rr')s(rr')T the product of singlet bonds between r and r' for both S and T, and by 

Qrr' = [P^' ~ l}[Prr> } ~ 1]. ™ find that H ST = -JE<rr>> Qrr' and 



4 



Q rr ,{rr')=A{rr') 
Q r n(rr'){ll') = (r'O(rZ') 



(14) 



First observation: the subspace of valence bond states in which bonds are simultaneously S- and T-singlets is left 
invariant by the hamiltonian. (In all figures, we will denote by a line such simultaneous S- and T-singlets.) Such 
a subspace is, however, still overcomplete. [|15| In D = 1, nevertheless, if we apply i?sT any number of times to 
the valence bond solid in Fig. [j], we generate a subset of valence bond configurations which are linearly independent 
(although not orthogonal). These are the valence bond configurations associated to non-crossing Lewis diagrams, fill 
It is known that such a basis can be conveniently used to work out efficient exact diagonalization algorithms in small 
chains. [[l6) The ground state wave function for a chain of N = 8 sites, for instance, can be easily obtained by solving 
a simple 3x3 problem. Fig. ^ shows the result obtained for the ground state wave function, and the lowest excited 
state of momentum 7r, in terms of VB configurations. 

Second observation: as is apparent from Eq. [l4|, a nearest-neighbor bond contributes a diagonal energy of —(4 J), 
whereas the off-diagonal matrix elements connecting different VB configurations are smaller by a factor 4, i.e., — J. 
For the SU(n)-mode\ of Affleck, the off-diagonal matrix element are smaller by a factor 1/n, which implies that only 
diagonal energies are retained in the limit n — > oo. The ground states, for n = oo, are therefore the VB configurations 
with the maximum number (= iV B ites/2) of nearest-neighbor bonds connected by VBs. In D = 1, there are only two 
such states, i.e., the VB-solids in Fig. nl 

To check if the VB scenario remains correct for our n — 4 case, we resort to exact Lanczos diagonalizations of chains 
up to 14 sites, and to Green Function Monte Carlo for longer chains. Fig. || shows the finite-size value of the gap 
between the ground state and the first excited state both in the Sj. OT = T^ OT = 1 (□) and in the S^ot ~ ^tot = 
sectors (CO- The full symbols refer to Hst, whereas the open symbols refer to two decoupled Heisenberg chains (no 
quartic term in Eq. ^|). All these excited states have momentum tt relative to the ground state. For the Heisenberg 
case, the excited state is a triplet and the model has gapless excitations with gaps decreasing as 1/L for L — > oo. 
jl7j For i?sT, the situation is different: The lowest excited state is a singlet, with a finite size gap AE(S = 0) 
decreasing faster than 1/L for L — » oo. The triplet, instead, lies above with a gap AE(S = 1) extrapolating to 
a finite value as L — > oo. This is a clear signal of a degenerate infinite volume ground state with a spontaneous 
breaking of translational symmetry and a gap to all excitations, consistent with a VB phase. Further support to this 
interpretation can be gained by direct inspection of the spin-spin correlations, (SfS^), the dimer-dimer correlations 
along a chain ((Sf S* +i )(S* , and on different chains ((Sf Sf +1 )(TjT 3 z +1 )) . Fig. || shows a log-log plot of the size 

scaling of the peak of the different structure factors, located at momentum tt. While the dimer-dimer structure factor 
peaks, ^(-^((SoSTXSfSf+x)) and Y, l {- l ) l {{ s o s i)( T ? T ?+i)), tend to diverge linearly with the length L of the 
chain, the usual spin-spin structure factor tends to a finite value in the thermodynamic limit. 

Finally, simple variational arguments also point to a VB ground state in D = 1, while being, we believe, not 
conclusive in D = 2. A VB phase benefits highly from the large values of (S r • S r >) = (T r • T r /) = —3/4 on the singlet 
bonds, due to the presence of the "square" of the singlet contribution in the quartic term of i?sT, which enters with 
a factor 4. Let us compare, for instance, the energy of the following two states: (i) The product of two independent 

Heisenberg ground states, |* 2 -h) = |*h ) ® l*H )' and (") tne product of two VB states, |* 2 -Vb) = I*vb) <H> I*vb) 
where each VB state has the simple form of a product of dimers on adjacent sites. In D = 1 we take 

|*vb> = (12) (34) ••■ (L-1L) , (15) 

where (ij) = (f J. — l^)/^/2 denotes a singlet between sites i and j. In D = 2 we consider the simplest short range 
VB state: a columnar dimer state, or, equivalently, any state with nearest-neighbor pairs coupled into singlets in an 
arbitrary manner (there is a huge degeneracy). The energy per site of the two states are: 

e 2 -H = 2e H - UD(e K /DJ) 2 , 
e 2 _VB = 2e V B - 4 J(l/2)(-3/4) 2 . (16) 

Here e^/J = D(S r ■ S r >) (r and r' nearest neighbors) and eyB = (— 3/8) J are the energy per site of the Heisenberg 
ground state and of the VB state, respectively. For the Heisenberg model we have en/J = (— In 2 + 1/4) w —0.4431 
(in D = 1) and en/ J =~ —0.66 (in D = 2), whereas for the crude VB state considered here, the expression for £2-vb 
turns out to be independent of D. The factor (1/2) in the expression for e 2 -VB is due to the fact that only half of 
the bonds in the VB state enjoy the large singlet-singlet pairing — 4J(3/4) 2 . In spite of this reduction, the coefficient 
of the quartic term is substantially larger for the VB state than for the Heisenberg state in both D = 1 and D = 2 
(« -0.2812(4J) for the VB state, to be compared to « -0.1964(4J) and « -0.218(4J) for the D = 1 and D = 2 
Heisenberg, respectively). In D — 1 the VB state wins over the Heisenberg state. Generally speaking, an increase 
in the coordination number tends to stabilize a Neel-likc antiferromagnet with respect to VB states. In D = 2, the 
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previous crude variational estimate would indeed give the Heisenberg state as favored. The VB state considered here 
is, however, very poor: for instance, its energy per site, neglecting the quartic term, is only — 0.375J, whereas it is 
well known that short-range RVB states close in energy to the Heisenberg ground state can be constructed. Q As a 
consequence, we believe the present variational estimate to be not conclusive in D = 2. 

In case the VB picture should be energetically favored, we face the problem of determining the possible orderings 
of the dimers in the lattice. In principle, either a VB crystal with broken translational symmetry or a homogeneous 
spin liquid may occur. Following Affleck's D = 1 approach, jl2) a possible way of tackling this problem, in D = 2, is 
from the n — > oo limit. For n = oo, once again, all the states made up of the maximum possible number (= iV s ;tes/2) 
of nearest-neighbor S- and T-VBs - we refer to them as NNVB configurations - are degenerate, linearly independent 
(although not orthogonal) ground states. Contrary to the D — 1 case, where only two such states exist (see Fig. [j]), 
in D = 2 this implies a huge degeneracy. For n — > oo, however, there are off-diagonal matrix elements connecting 
different NNVB states (smaller by a factor 1/n, as previously discussed), which lift this degeneracy to first order. It 



is straightforward to verify, using Eq. 14, that the effective hamiltonian within the subspace of NNVB states (to first 



order in 1/n) is just the quantum dimer model (QDM) first introduced in Ref. |19| 

Hqmd = — J {I =>(lll+H.c.} + y Yl {I =>(= l + IIIXIII}. ( 17 ) 

plaquette plaquctte 

for the particular case of V = 0. (Here || and = schematically indicate adjacent S- and T-VBs in the y and x direction.) 



The QDM is known to have a columnar solid phase in a region of parameter space close to V = 0. |19 20] Recently, 
an exact diagonalization study of the QDM has revealed that such a columnar solid phase has, more precisely, the 
features of a plaquette RVB, i.e. translational invariance in one direction is spontaneously broken, but the spin-Peierls 
"columns" are disordered as a consequence of a resonance of parallel bonds in each plaquette. pi| This picture would 
yield a fourfold degenerate ground state in two dimensions. 

We conclude that, for n — > oo, our model has a columnar VB phase in D = 2. Given the fact that our crude 
variational calculation appears to favor a Neel state in D = 2, we believe that the issue of the actual behavior of our 



n = 4-case, H$t, in D — 2 is still open. This issue will be addressed in a future numerical study. 22 



IV. LESS THAN HALF-FILLING. 



The motion of holes in a VB state has been the subject of considerable research starting from the early works on 
high temperature superconductivity J23|. Now, we are able to investigate this problem on the basis of a consistent 
microscopic model which does have a VB solid as a ground state. This study may shed light on the issue whether 
holons and spinons decouple in VB states and whether the VB spin background favors effective hole-hole attraction. 
Therefore, we analyze in some detail the case of filling n < 1 in one dimension. For large repulsive interactions (small 
J) the effective hamiltonian describing the system is simply the analog of the t — J model 

H t _ ST = P KP - 4J ^ n r n r >{S r ■ S r < - ^)(~T r • T r < + 2T ? ;T*, - -) + (3 - sitctcrms) , (18) 

<rr'> 

where K is the hopping part of the original hamiltonian (Eq. |]), P a is the projector onto states without double 
occupancy, and n r is the occupation number operator. In D = 1, analogously to the Hubbard case, p5|p^ | the 
(real-space) wavefunction factorizes, for J — > 0, into a spinless fermion part times the wavefunction of a suitable spin 
hamiltonian H spm on the squeezed chain |27]] 

• • • ,xjv;cri, • • • ,erjv,Ti ■'-,tn)= iPsf(xi, ■ ■ -,x N )$(o-i, ■ ■ -,a N ,Ti ■ ■ ■ ,t n ) . (19) 

The position of the electrons is determined by the spinless fermion wavefunction ^SF^i, • • • , xn), while the spin and 
orbital ordering is governed by a spin hamiltonian which is obtained by taking the average of H t -ST over the spinless 
fermion state. The result is simply related, by an overall change of energy scale and a canonical transformation, to 
the previous spin model Hst (Eq. pi] ): 

ffspin = -4Jn[l - SJ ^^} E ( S * ' S ^ - i)(- T * ' T -> + 2T * T ! ~ \) ■ ( 2 °) 

<ij> 

The physical origin of this factorization relics on a general feature of one dimensional models with nearest neighbor 
hopping: In the U — > oo limit, the spatial ordering of the particles is conserved by the dynamical evolution, irrespec- 
tively of the spin configuration. Therefore, the hamiltonian can be diagonalized in each subspace defined by a given 
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spin ordering and the result does not depend on the chosen spin configuration. The effective hamiltonian governing 
the charge degrees of freedom is that of a free spinless fermion gas (with suitable boundary conditions which do depend 
on the spin wavefunction psfl) whose eigenfunctions are ipSF- The degeneracy among different spin configurations is 
lifted to O(J) by the magnetic term in H t -sT (0), leading to the complete classification (|l9| ) of the eigenstates of 
H t -sT to lowest order in J. 

The spin hamiltonian in Eq. |2^, we argued before, has a spin gap. Therefore, we conclude that, while the charge 
sector is gapless (metallic) for large positive AEs and n / 1, the spin sector is gapped at all densities and has a 
doubly degenerate ground state in the thermodynamic limit. However, such a spin gap is not related to the formation 
of bound electron pairs. In fact, the charge degrees of freedom behave as free particles with no tendency towards 
pairing: the superconductive correlations are not enhanced in the ground state and decay as ~ x~ 2 . This state can 
be rather interpreted as the one dimensional realization of an "itinerant" valence bond state where spinons bind, 
leaving holons in a Luttingcr Liquid state. (The situation is quite similar to the strong coupling phase of the t — t' 
Hubbard model in D = 1, for t' > t/2. j24|]) The very presence of such a phase is then a remarkable consequence of 
one dimensional spin-charge decoupling. 

The factorization of the wavefunction at J —> can be used to calculate various physical quantities of this spin- 



gapped metallic phase, following analogous calculations in the Hubbard case. [|25 28 1 In particular, density correlation 
functions just coincide with those of a spinless fermion gas, showing power law decay at large distances (~ x~ 2 ) 
while the spin-spin (and pseudospin) correlation still decays exponentially, similarly to the half filled case. Another 
interesting property which can be evaluated is the one-hole Green's function providing a quantitative description of the 
motion of one hole in a short range VB state. The picture of a static VB crystal would yield an exponentially decaying 
Green's function in real space, but ID fluctuations severely affect the asymptotic form of the Green's function. The 
calculation can be carried out exactly in the J — > limit, where the factorization property of all eigenfunctions holds. 
In fact, in this limit, the energy levels of a single hole in a VB state coincide with those of a free particle in vacuum 
plus O(J) corrections depending on the energy of the spin background. Moreover, the energy, to lowest order, just 
depends on the holon momentum, i.e. on the momentum k of the spinless fermion part of the wavefunction ^sf- 
E(k) = 2tcosk + O(J). On the other hand, the holon momentum is related to the total momentum p of the state by 
the conservation law: p = k + Q. Here Q is the momentum of the spin part of the wavefunction $ and is quantized 
in units of 2tt/ (L — 1), the spin model being defined on the squeezed chain of L — 1 sites. Therefore, for J — ► 0, the 
one-hole energy levels of the model, at a given total momentum p, can be written in terms of the spinon momentum: 
E p Q = 2 cos(p— Q), i.e., of the momentum of the spin part of the factorized wavefunction. The corresponding density 
of states coincides with that of a single free fermion: 

N ^'^7WTW' (21) 

with singularities at the band edges to = ±2i. These singularities are due to the quadratic dependence of the energy 
spectrum on the holon momentum which will be probably cut-off by the presence of spin excitations to O(J). 
The calculation of the single hole Green's function 

G( P , t) = i (*| C t )CT e -it{H-E -iS) Cp a W 0{t) (22) 

is conveniently carried out by use of the Lchmann representation of the corresponding spectral function: 

A(p, u ) = ilm G(p, u>) = V | (tf * | c v , a | *) | 2 5{w + E - E.) 

7T * — ' 

s 

= Y,Z{Q)K"-E p ,q), (23) 
Q 

where, in the last equation, use has been made of the special form of the energy spectrum as J — > 0. The usefulness 
of this representation rests on the knowledge of the weight function Z(Q) which, following Rcf. p8| |, can be written 
in terms of a non-local correlation function involving only spin (and pseudospin) variables: 

L-2 



z(Q) = \ j^Y. e ~ lQ3 ^) 



3=0 

= ($ | (2Tj ■ Tj-i + i)(2Sj • S^i + i) . . . (2Tj • T + i)(2S x ■ S + ~) | <&> . (24) 

Such a correlation function can be evaluated on the ground state of the H^t hamiltonian by Lanczos diagonalization 
and, similarly to the Heisenberg model case, shows negligible size dependence. Calculations have been carried out up 
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to L = 18 sites, and the results are shown in Fig. 0. A clear singularity in the quasiparticle weight Z(Q) is visible at 
the spinon Fermi momentum Qf = tt/2. The analysis of Lanczos data suggests a power law divergence as Q — > Qf> m 
close similarity to the analogous behavior encountered in the gapless Heisenberg model. A quantitative determination 
of the critical exponent is however precluded by the severe limitation in system size associated to the Lanczos method. 
In the thermodynamic limit, this singularity gives rise to branch cuts in the analytic structure of the spectral function 

A(p,lo) = [L - 1) [Z (Q + (p,u)) + Z(Q-(p, u>))] N(u>) (25) 

along the lines defined by \Q±\ — n/2, where Q±(p,u>) — p ± arccos(^) i.e. at ui p = ±2sinp. The form of this 
dispersion relation suggests nearest neighbor hopping processes of the hole across the lattice. 

The overall picture of hole motion in a Valence Bond state is therefore quite similar to the case of holes in one 
dimensional quantum antifcrromagnets, being characterized by power law tails in the single particle Green's function 
at the dispersion energy uo p . This behavior can be traced back to the presence of branches of gapless collective modes 
above the spin gap which allow the decay of the spin 1/2 excitation created by the hole motion. These features of 
hole propagation should be contrasted to the properties of charges in Ising antiferromagnets, where the excitation 
spectrum is discrete and the band is dispersionlcss. 

V. CONCLUSIONS 

In this paper we have discussed in some detail a hopping model of strongly repelling electrons in the presence of 
local orbital degeneracies not removed by crystal field effects. 

For the insulating case (one electron per site and large Coulomb repulsions) we have studied a particular exchange 
hamiltonian, Eq. (fil) . It describes the low-energy physics of this system in the regime in which the inter-orbital singlet 
is the lowest-energy intermediate state available to hopping. This regime tends to be further favored by a relatively 
strong electron-phonon coupling. At variance with the standard Hund's rule case (the high spin two-particle state 
has lowest energy), which would instead favor spin ferromagnetism, we find in our case a strong tendency to the 
formation of a statically ordered Valence Bond (VB) phase. In D = 1, this has been firmly established by using 
variational arguments, a large-n limit, and exact diagonalization results. Given the presence of a sizeable spin-gap, 
weakly coupled chains should also have a VB ground state. 

More uncertain is, instead, the outcome in the isotropic D = 2 case, where a large-n approach would predict a 
columnar VB phase, whereas a crude variational calculation rather points to a Neel-like antiferromagnet. A numerical 
study, presently under way, will clarify, we hope, this matter. 

Finally, we have considered the behavior of the system upon hole doping in the D — 1 case. The factorization of 
the wavef unctions, valid in the strong repulsive limit and peculiar to the D = 1 case, is the crucial ingredient used to 
tackle the problem. The outcome is a kind of "itinerant" VB state, where charge degrees of freedom are gapless, but 
the spin and orbital degrees of freedom are gapped and VB-like. 

The motivation for our study has been primarily theoretical. Originally inspired by dynamical Jahn- Teller molecular 
lattices, such as Cqq compounds, the ingredients used are, admittedly, extremely idealized, and not easily applicable 
to an existing realistic situation. In a real JT lattice, both the crystal field effect, and the intermolecular elastic 
coupling cannot normally be neglected, and may favor a cooperative static JT distortion of all molecules. Q Second, 
intermolecular electron hopping tends to be comparable, or even larger in comparison with relevant vibron energies. 
Apart from this, realistic hopping matrix elements have a strong orientational dependence in more than one dimension, 
and, in general, there is no reason for neglecting overlap integrals between orbitals with different orbital quantum 
numbers (see Eq. (|4|)). Third, a detailed description of the structure of the doubly occupied states at each site brings 
in correction terms to the exchange hamiltonian which turns out to be, in general, strongly asymmetric in orbital 
space (see Eq. (|To|)). 

On the opposite side, the exchange hamiltonian we have studied is elegant, parameter-free and prototypical of 
strong correlations in an orbitally degenerate lattice. The physics uncovered has a definite robustness against "small" 
perturbations, due to the gap present in the excitation spectrum. The VB gapped state has a clear resemblance to a 
spin-Peierls state in D = 1. The spin-gapped metal found upon doping is interesting. The basic ingredients needed 
for a realistic system to be a candidate VB described in this work are: i) orbital degeneracy not trivially removed by 
cooperative JT and/or crystal field effects; ii) relatively large molecules with a strong electron-phonon coupling, so 
as to make the inter-orbital singlet favored as compared to the Hund's rule triplet; iii) narrow bands with relatively 
large on-site Hubbard U, so as to stabilize a Mott-Hubbard insulator at exactly one electron per site; iv) reduced 
dimensionality. 
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Ceo charge transfer compounds, based on Cgg ions, EJ^] are potential candidates to the realization of such a 
scenario. Cqq has a triply degenerate t\ u molecular orbital, coupled to several intramolecular vibrations, resulting in 
an important dynamical JT effect. The electron-phonon coupling leads to a substantial pairing energy w O.leV" 
which is, however, overwon by a substantially larger Hubbard U sa 1 — l.SeV^, [ j30[ resulting in a Mott insulating 
behavior. The lattices tend to have anisotropic lattice constants, with a pronounced quasi-one-dimensional character. 
p9f As for the virtual intermediate states of Cqq , singlets and triplet are close in energy, at least in solutions. (3^] If 
the triplet prevails, as is possibly the case in TDAE-Ceo, f|,|| there can be spin ferromagnetism. We believe, however, 
there is room in some future system for the alternative possibility of a spin gap, and the formation of the spin-orbital 
VB state described here. 
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FIGURE CAPTIONS 



FIG. 1. The two VB solid configurations in D — 1. 



FIG. 2. The ground state (momentum 0, +) and first excited state (momentum n, — ) for Hst on a chain of 8 sites in terms 
of S- and T-VB configurations, (a — 0.566 • • •, j3 = 0.416 • ■ •) and (a = 0.268 •••,/? = 0.102 • • •), respectively, for the two states. 

FIG. 3. Finite size gaps for the lowest excited states of the hamiltonian Hst in D = 1 (solid symbols) compared to the 
Heisenberg chain case (open squares). The data are obtained by exact diagonalization, for chains up to 14 sites, and by Green 
Function Monte Carlo for longer chains. The excited state for the Heisenberg case is a triplet and AE/J = 7r 2 /(2L) + ■ ■ ■, 
shown by the dashed line. The solid straight line, extrapolating to a finite value for L — » oo (~ 0.5969), is obtained by a 
quadratic least-square fit to the triplet excitation results. This figure shows that, whereas the triplet excitations are gapped, 
the singlet gap decreases faster than 1/L, signalling spontaneous symmetry breaking, consistent with a VB ground state. 



FIG. 4. Log-log plot of the 7r-component of the Fourier transform of different correlation functions for the hamiltonian Hst in 
D — 1. Open squares refer to (SfSj), open and solid circles refer to {(Sf S', 2 +1 )(S|S| +1 )) and {(S* Sf +1 )(T/T/ +1 )), respectively. 
The dashed line has slope 1, for comparison. The dimer-dimer correlations grow fastest with size, again consistent with a VB 
ground state. 

FIG. 5. The function Z(Q) defined, in Eq. 0, determining the single-hole spectral function in the thermodynamic limit (see 
Eq. E5J). The data are obtained from exact diagonalizations of chains up to L — 18. 
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